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Abstract 

We investigate the conditions under which systems of two differential eigenvalue equa- 
tions are quasi exactly solvable. These systems reveal a rich set of algebraic structures. 
Some of them are explicitely described. An exemple of quasi exactly system is studied 
which provides a direct counterpart of the Lame equation. 
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1. Introduction 

Quasi exactly solvable (QES) equations are differential, spectral equations (typi- 
cally Schrodinger eigenvalue equations) for which a finite number of eigenvalues and 
eigenvectors can be found by algebraic methods. In this respect, they contrast with 
the so called completely solvable equations (like the quantum harmonic oscillator or the 
Coulomb problem) where the whole discrete spectrum can be computed algebraically. 

After their introduction, a few years ago, by A. V. Turbincr [1] and A. G. Ushvcridzc 
[2], the quasi exactly solvable equations received many new developments (see Refs. [3,4] 
for reviews). 

Scalar QES equations in one variable are now well understood [1-5]. The basic 
ingredient is that, in a good variable, the linear differential operator defining them 
can be expressed in terms of a projectivized representation of the generators of the 
group SL(2,IR ). Scalar equations in two variables, one real and one Grassmann, were 
investigated in Refs. [3,4,6]. All results obtained in this context can be reformulated in 
terms of systems of two equations in one real variable. The relevant algebraic structure 
behind this type of equations is the graded algebra osp(2,2). 

In this paper, we want to revisit the systems of two coupled QES equations from a 
more general point of view. The understanding of these equations requires the knowledge 
of all 2x2 matrix differential operators that leave invariant the doublets of polynomials 
of given degrees, say m and n. Such an investigation is reported in Sect. 2, where we 
put the main emphasis on the algebraic aspects of these sets of differential operators. 

In Sect. 3 we describe the possible forms of systems of two quasi exactly and com- 
pletely solvable equations. We compare the number of free parameters they can depend 
on. Finally, in Sect. 4 we present an explicit example of such system. It is related to the 
stability of sphaleron solutions in the two dimensionnal Abelian-Higgs model. 

Before to start with the systems, we briefly recall the basic features of the scalar 
QES equations in one real variable. They are strongly related to the linear, differential 
operators that preserve Vn- the space of polynomials of degree n in a variable, say x. 
Among these operators, the three given by 



play a considerable role. They obey the commutation relations of the generators of the 
group SL(2,IR ): 




(1) 




[J, 




(2) 
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and therefore constitute an irreducible representation of dimension n + 1 of the algebra 
of this group. Using the irreductibility of the representation it is easy to show that the 
linear differential operators preserving the space P„ Sire the elements of the envelopping 
algebra generated by the three operators (1). Accordingly, if P denotes a polynomial 
in three variables, then the eigenvalue equation 

P{J+,4,J-)-F{x)^XF{x) (3) 

is equivalent to a system of n + 1 algebraic equations. If P is quadratic, it is possible, 
formally, to rephrase Eq. (3) as a Schrodinger equation by mean of a change of variable 
and of a change of function. 

2. Algebraic Structure of QES Systems 

The purpose of this section is to describe the algebraic features of quasi exactly 
solvable systems of two equations. We denote by Vm,n the space of doublets of polyno- 
mials in X, the first (resp. second) component being of degree m (resp.n). Given the 
integers m, n, we will consider the set of operators defined by 

rp+ ^ ( Jm ) rpO _ ( Jm 1 rp- _ ( J-m 

^ " ' J+ j ' ^ - I JO J ' ^ - 1^ J - 



, 1 fn + A , . _ 
J = - „ ,A = n — m 



2V n, 

n.-ce , , 

Q^ = a;"-V- , a = l,2,---A + l 

It is easy to show that these operators are linearly independent and that each of them 
preserves the space Vm,n- Moreover, reasoning component by component, one can 
convince that any linear operator preserving Vm,n can be constucted polynomialy from 
the 6+2 A generators (4). So, the operators (4) play for systems of two equations the 
same role as the operators (1) do for scalar equations 

We want to demonstrate that the operators (4) close under (anti-) commutation, 
independently of m and n. The subalgebra sustended by T^,T^ is isomorphic to the 
algebra of SL(2,IR). The SU(2) generators, say Tf. can be recovered by the combinations 



Ti = i(r--T+) , T2='-{T-+T+) , T3 = T° 



(5) 
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The commutation relations and the Casimir operator associated to this (reducible) 
representation are respectively 

[Tj,Tk] = ie-ikiTi (6) 



The operators (4) behave differently under the commutation with the matrix J 
(whose form in (4) is defined for later convenience): 

[J,Tfe]=0 , [J,Qa] = -^Qa , [J,Qa] = ^Qa (8) 

From now on we will refer to T^, J as to bosonic operators and to Q (resp. Q) as to 
fermionic (resp. anti-fermionic) operators. The J-weight of bosonic operators is zero 
while that of fermionic (resp. antifermionic) ones is —A/2 (resp.+A/2). 

For the commutation rules between bosonic and fermionic operators we obtain 

[Qa,T+] = (l-a + A)g«+i 

[Q«,tO] = (1-«+^)Q« (9a) 
[Q«,T-] = (l-a)g«-i 

[Qa,T+] = -{l-a)Qa-i 

[Q«,T«] =-(!-«+ |^)Q« (96) 
[g^T-] = -(!-« + A)g«+i 

demonstrating that the sets of Q^s and Q^s transform according to the s = ^ repre- 
sentation of the SL(2,IR ) subalgebra. This statement becomes more transparent if we 
redefine these operators according to 

2s V^' . . A 



The computation of the anti-commutation relations between fermionic and anti- 
fermionic operators is more involved. Formally, we obtain 

{Qa,Qp} = Map{T-r-'' ,iia>p 
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with 

A-a 13-2 

Map=\[{T'' + J,-j^-{l3-\)P2) \[{T'' + J-k^-{A + \-a)P^) (12) 
and 

Je=(A-l) I-J, Pl= l\,P2={^, ?l (13) 



^0 oy ' ^ vol. 

In the definitions (12) it is understood that the first (resp. second) product is I if 
O! > A (resp. if /3 < 2). The form (11) is not satisfactory because the projectors 
Pi, P2 entering in M.ctj3 are not expressible in terms of J and of C (the Casimir) in an 
n-independent way. These matrices however obey the identity 

(J + (^^))P2 = -^iJ" + i^- 2A) J + A(A -l)-C) (14) 
In order to demonstrate that Aiap is independent on n, Pi, P2 we write formally 

Map = AapP2 + Pa/3 (15) 

Using the fact that Pi, P2 are projectors, a direct evaluation of the matrix A^p yields 
Aap = \[{T' + X-j) n {T'-X-3)-\[{To-X-j) n (To + X-j) (16) 
with 

X^J+{^^) (17) 

and . . . . 

Q= 2 'g = P~Q'+ 2 >oi = -Qi+ 2 '"2 = P ^ (18) 

From Eq. (16), it is clear that the matrix A contains X as a factor. Then the use of (16) 
together with the identity (14) demonstrates the statement that M. are independent on 
n and on the P's. For completeness we also mention the relations 

{Q«,Q^} = ,{Qa, Q/3} = (19) 



So, Eqs.(9)-(19) finally show that the 6+2A generators (4) close under (anti)- 
commutation to a family of operator algebras indexed by the integer A. In order to find 
the abstract algebras for which the above operator algebras provide the representations 
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we have to check whether all Jacobi identities are fuUfiUed on abstract level and, if not, 
to modify some of the commutation rules. 

It appears that only the Jacobi identities for Q,Q and Q or for QQ and Q impose 
further constraints on generators. This suggests that just the anticommutation rules 
(19) should be changed. The anticommutator {Qa, Q/?} has weight 2 with respect to J. 
Therefore its value should also be quadratic in Q's (we exclude the more complicated 
objects like QQQQ etc.). By analysing the Jacobi identity for QcnQpiQ-y we arrive 
at the conclusion that the coefficients of the quadratic form on the right hand side 
cannot depend on T^'s and on J; i.e. they have to be numerical. Therefore the only 
possibility is that we decompose {Qa, Q13} into SU(2) irreducible components and keep 
only some of them on the right hand side. We solved this problem for A = 2 and 3. 
It appeared that in order to satisfy the Jacobi identity it is sufficient to keep only the 
lowest non vanishing representation. The resulting relation can be viewed as a quadratic 
constraint on the Q's. After this constraint is imposed, no further constraint is implied 
by the Jacobi identities. 

We now collect the cases where an abstract algebra was found. 

Case A = 

In the case A = 0, we have an ambiguity in defining the algebra. We can treat it 
according to the general scheme above (in this case J oc I). Alternatively, we can also 
put J = a"3 and impose the commutation rule on Q, Q. In this case one obtains the the 
Lie algebra of SL(2,IR )xSL(2,IR ) 

Case A = 1 

This case was considered at lenght in Refs. [3,4,6]. The generators (4) constitute a 
representation of the graded algebra osp(2,2). The anticommutators of Q and Q close 
into linear combinations of and J (see Eq.(9)). So, the relevant algebra is still a Lie 
(super) algebra; this is not the case anylonger when A > 1 



Case A = 2 

The algebra associated with the case A = 
operators in the vector basis, i.e. 



2 is better handled when we express the 
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^ Q3 - Qi ^ Q3 + Qi ^ ^ 
^1 = — K — ) '^2 = — :r. — , = -Q2 

- - _ _ (20) 

^ Q3 - Ql r\ Q3 + Ql r\ r\ 
C'l = ^ 5 C'2 — , C'3 — V2 

In this basis, the relations (8), (9), (10) become respectively 

[J,G,] = -G, , [J,e,]=G, (21) 

[Ti, e,] = zeijfcefe , [Ti,e,] =ieyfc0fe (22) 

{e„ e,} = 5,,(^ + ^^y^) - ^(T,T, + T,T, - y>,,C) + ie,,fe(J - \)T^ (23) 

In the last relation, the right hand side is decomposed according to irreducible represen- 
tations of SU(2). Using such a decomposition for the anticommutators (19) and keeping 
only the scalar part on the right hand side, we obtain instead 

{e„ea = ^5,,(e? + ei + ei) (24) 

and a similar relation for G. We checked that all the Jacobi identities are fuUfiUed with 
the choice (24). 

Case A = 3 

In this case, we found it convenient to express the operators Qa, Qa = 1, 2, 3, 4) in a 
spin-tensor basis. It consists in linear combinations of the Q, say Qm, where the index 
i = 1, 2, 3 (resp. a=l,2) transform according to the s — \ (resp. s — 1/2) representation 
of SU(2). The components of the spin-tensor Qj^ read 

Qgi = V2g3 , Qgi = ^2(53 
Q3I = V2Q2 , 03-1 = -^Qi 



Qii =i=(Q2-Qi) , Qi-^ = ^(Q2-Q4) 
Qi_i =^(gi-Q3) , Qi_i =-^(Qi-Q3) 



(25) 



g2-l = ^(Ql + Qs) , Q2-i = ^(<5l-Q3) 



and obey the identities 

Qia{cri)ab = , {(7i)abQib = (26) 

which guarantee that only the components corresponding to the representation s = 3/2 
are selected. 

The advantage of the spin tensor basis is that the commutators (9) write in a 
manifestly covariant way: 



(27) 



as well as the anticommutators (11): 



2 

{Qia, Qjb] = 2Tj(r • a)abTj - -C6ij{T ■ a)ab 

i 1 
+ -^^ijkTkCSab — -^C{Tiaj + Tjai)ab 

+ 2{J-l){TiTj-^5ijC)5ab 
3 

+ ~ 2^^'^jk{Tk, Ti}{ai)ab 

+ 2^{'^j^'^k}^ik£ — {Ti,Tk}ejk£)i(T£)ab (28) 

+ ^iJ-2){5J-3)eijkTkSab 



1 

3* 

1 



+ - J(J - 2){Ti{aj)ab + Tj{ai)ab - 45^ (T • a)ab) 



+ -J{J ~ 1)(J - 2){itijk{(Jk)ab - 26ij6ab) 

Keeping again the lowest dimensional representation appearing in the decomposi- 
tion of anticommutator {Qia^Qjb} (in this case it is the vector representation because 
the scalar piece identically vanishes), i.e. 

1 2 

{Qia,Qjb} =T:io'2)cd{QidjQjc}io'2)ab + -^5ij{Qka, Qkb} 

^ ^ (29) 



-^{{P'2(yi)ab^3kl + {(^7.0 j)ab^%kl^{p-2)cd{(ikd-,Qkc) 
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we have checked that all Jacobi identities are fuUfiUed. 

3. Systems of QES equations 

The 2x2 differential matrix operators which are quasi exactly solvable can be ex- 
pressed in terms of elements of the envelopping algebras constructed from the generators 
(4) . Using the same notations as in Ref [6] , these operators have the following form 

Tfe(x) = J]afe,,(x)— (30) 

where the coefiicient ak,i{x) are 2x2 matrices with polynomial entries. Their degrees 
are indicated here in the square brackets 

/ ^[k+i-A]\ 

(^k,i{^) = ^[k+i+A] j^[k+i] I (31) 
\'-^k,i ^k,i I 

Off course not all the polynomials Ak,i,Bk,i ■ ■ ■ are arbitrary: for k > A and for 
generic n, the most general operator Tk{x) depends on 

4ik+lf , iik+lf) (32) 

arbitrary parameters, independly on A (the corresponding number for scalar equations 
is given in parenthesis). For k < A the situation is more complicated because some 
constraints may become dependant. 

For the operator Tk{x) to be exactly solvable, the degree of the polynomials in ak,i 
need to be as follows 

-M(-) = L[i';] (31') 
\^k,i ^k,i / 

they are not constrained, so that these equations depend on 

> A 

(32') 



2{k+l)ik + 2) + ^^ — ^ for k>A 



^{k+l){k + 2) + {k+l)A for k<A 



free parameters (this number (A; + 1)(A; + 2)/2 for scalar equations). 
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4. Example 

In this section, we discuss a system of two coupled equations which admits algebraic 
solutions. This example arises in the study of the stability about sphalerons [7] (i.e. 
unstable classical solutions) in the Abelian gauge-Higgs model in 1+1 dimension [8]. 
The relevant Schrodinger equation reads 



^ + A - e^kW -2dk cn • dn \ / f{z) 

-2dk cn • dn ^ + A + 1 + /c^ - (6>2 + 2)kW J \ 9(4 



(33) 



and is considered on the Hilbert space of periodic functions over [0,4K(k)] (K(k) is 
the complete elliptic integral of the second type). The three elliptic functions [9] 
sn{z,k),cn{z,k),dn{z,k) are periodic with periods 4K(k), 4K(k), 2K(k) respectively. 
The spectral parameter A represents the mode eigenvalue: negative A's correspond to 
the unstable modes of the sphalerons. 

The parameter depends on the coupling constants of the model; it represents the 
mass ratio 2Mh/Mw where Mw (resp. Mh) is the mass of the gauge (resp. Higgs) 
boson. The system above admits algebraic solutions if these masses are such that 

= N{N+1) or M| = ^^^^ "^W ^ with N integer (34) 
Following a similar scenario as for the Lame equation [10] 

(^ + A - N{N + l)khii'')f{z) = (35) 

the algebraic solutions of Eq.(33) occurs in four sectors of the Hilbert space [11]. In 
order to construct them, we need to set the system in the forms discussed previously. 
For this purpose, we perform the change of variable 

X = sn{z, kf (36) 

accompanied by a change to new functions P{x), Q{x) defined through 

f{z)^F{z)P{x) , g{z) ^ G{z)Q{x) (37) 

where F{z)^G{z) are some products of sn, cn, dn. Then we have to determine F{z) 
and G{z) so that the new equations for P{x),Q{x) admit polynomial solutions in the 
variable x. The factors F{z)^ G{z) allowing for an algebraisation are different according 
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to the parity of N; the eight different possibihties are reported in the first two columns 
of Table 1. 

Even if the changes of variable (36) and of function (37) are performed, the equation 
for P, Q is still not in a form that preserves Vm,n'- an additional change of basis of the 
form 

need to be done. In the eight cases, the change of basis V can be constructed from the 
following matrices: 

They are presented in the third column of Table 1 and the dimensions m, n of the 
invariant space Vm,n are reported in the last two columns, these results demonstrate 
that Eq.(33) admits four types of algebraisation, yielding a total of 4N + 2 algebraic 
solutions, if N is an integer (this number is 2N + 1 for the Lame equation). The first 
(resp. last) four lines in the table correspond to the solutions available for N odd (resp. 
even). 

The analogy between Eqs. (33) and (35) is also present when doubly periodic 
solutions [10] are considered, i.e. solutions of period 8K(k). In order to discuss this 
issue we consider the following changes of functions: 



f{z) = Vdn(z) ± cn(^)(±y(a;)cn(^) + Z{x)dn{z)) 
g(z) = y/dn{z) ± cn(^)(±F(a;)cn(^) + W{x)dn{z)) 



(40) 



and 



(41) 



f{z) = Vdn(^) ± cii{z){±Y{x)cii{z)dii{z) + Z{x)) 
g{z) = Vdn(z) ± cn{z){±V{x)cn{z)dn{z) + W{x)) 
Due to the square root factor, these functions have 8K(k) as period, the functions with 
+ and — correspond to each other by a translation z ^ z + 2K{k). 



Inserting the ansatzes (40) into Eq.(33) and identifying to zero the coefficients of 
the factors cn, dn, yields a system of four equations in F, Z, V, W. To put this system 
into a cannonical form; we must perfom an additional change of basis. In this case, it 
reads 

\ /^\ 

; (42) 





/ 


w 




z 






V 



1 

ke 

N+l 






1 

-ke 

N 









1 

ke 

N+l 






1 

-ke 

N 



w 
z 
\v/ 
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then the differential operator associated with the ansatz (40) preserves the space (using 
obvious notations) 

^iv-i,iv,iv-i,iv if N=^+2N, iV integer (43) 

Similarly, after the change of basis defined by 

/Y\ / 1 \ /Y\ 
W \ _ 1 I W 

Z - 1 \Z 

\v / V 1 / \v/ 

we obtain the operator associated with the ansatz (41) in a form that preserves the 
space 

^iv,iv,iv,iv if N=^+2N, N integer (45) 

This demonstrates that 4A'" + 2 algebraic solutions are available also when N is a semi- 
integer once the degeneracy ± in Eqs. (40,41) is taken into account (both signs yield 
solutions of equal eigenvalues). 

5. Conclusions 

Quasi-exactly-solvable equations constitute an attactive bridge between group the- 
ory and spectral equations. While ordinary QES equations are related to usual Lie 
algebra, systems of coupled QES equations are related to more sophisticated structures. 
For systems of two coupled equations, the relevant operators are those that preserve the 
spaces Vn,m- These operators can be perceived as projectivised representations of some 
graded algebra. For our equations, the relevant graded algebra have finite numbers of 
generators and our results demonstrate that they admit representations of arbitrary 
finite dimension. The anti-commutators between the fermionic generators close into 
n — m-powers of the bosonic generators. 

Applications of quasi exactly solvable systems can be found in the framework of 
quantum mechanics, namely in the study of coupled channels. Here, we illustrate the 
relevance of such systems within another domain, namely in the study of unstable modes 
about sphalerons in a simple field theory. This example suggests that the stability 
analysis about other classical solutions (solitons, kinks....) might also be related to 
QES systems (see, however, Ref.[12]). 



(44) 
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TABLE 1 



F(z) 


G(z) 


V 


m 


n 


1 


cn dn 


a B(l/k) 1 


N-l 
2 




N-l 
2 


cn dn 


1 


a (3{k) 7* 


N-l 
2 




N-l 
2 


sn cn 


sn dn 


a p{l/k) 6 


N-3 
2 


N-l 
2 




sn dn 


cn cn 


a I3{k) 6 


Af-3 
2 


AT 1 

2 




sn 


sn cn dn 


a l3{\/k) 7 


JV-2 
2 


7V-2 
2 




sn cn dn 


sn 


a P{k) 7* 


N-2 
2 




N-2 
2 


cn 


dn 


a (3{l/k) S 


N-2 
2 




N 
2 


dn 


cn 


a P(k) S 


N-2 
2 




N 
2 



The factors F{z),G{z) allowing for an algebraisation of the system (33) are given in the 
first two colums. The third column contains the change of basis (38) and the last one 
the degree of the invariant polynomial space Vm,n 



